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This paper provided a comparative study between linear and nonlinear 
controllers of a solar photovoltaic (PV) water pumping system using an 
induction motor and a centrifugal pump. For linear controller, classical 
Perturb and Observe was selected to ensure the operation of the PV system at 
the maximum power point (MPP) and is combined with Indirect Rotor Field 
Oriented Control (IRFOC) based on a conventional proportional integral 
speed regulator chosen to control and to optimize the rotor speed. In second 


Keywords: and third controllers, Backstepping and first order sliding mode controls were 
Backstepping control proposed for controlling the whole system. To regulate and optimize rotor 

i speed in the nonlinear controller, all considered control techniques were 
Centrifugal pump combined with IRFOC in order to establish easy control laws. In addition, 
Induction motor MPP was tracked by acting on the DC-DC converter and using its 
IRFOC mathematical model for developing control laws. Maximum delivered power 
Lyapunov functions was used as reference signal for optimizing actual rotor speed. The controlled 
MPPT system is operated without mechanical sensors. Estimators of rotor speed and 
P&O load torque were proposed based on the mathematical model of induction 
PV system motor and centrifugal pump and using only available output measurements.. 


Simulation results were investigated and the effectiveness of the nonlinear 


Slinding mode control : 
proposed strategies. 


Water pumping system 
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1. INTRODUCTION 

Solar photovoltaic water pumping system is a promising solution to solve the water shortage 
problem in Morocco. This is especially interesting in remotes areas, where the country has a considerable 
solar potential over 3000 hours of yearly sunshine and an incident solar radiation greater than 5 KWh/m? [1]. 
Furthermore, to meet its increasingly demand and to decrease its fossil fuel energy sources dependence, 
Morocco launched the Moroccan Solar Plan as part of its national energy strategy 2030. Due to this strategy, 
35 percent of its electricity from renewable energy sources had been produced in 2018 and it aims to produce 
52 percent by 2030. In the next ten years, the Kingdom has a plan to massively exploit this clean 
and intermittent energy and to reduce its energy dependence on foreign countries, from 95% to 85% by 2020. 
Solar photovoltaic energy offers an advantageous and competitive alternative compared to other traditional 
energy resources. According to Moroccan Agency for Energy Efficiency, these systems remain the cheapest 
energy sources and diesel is the most expensive energy sources [2]. Today, characteristics of economic 
growth based on the consumption of oil are measured in terms of ethics. Regarding this non-renewable fossil 
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resource, it took several hundred million years to produce it and only a few decades to exhaust it. If we add 
fossil fuel depletion and climate change, then we measure how much the energy transition is necessary and 
even urgent. The development of renewable sources is essential, with an associated absolute priority. The use 
of conventional oil energy in pumping system can disturb the livestock, pollute the agriculture environment 
and require transportation that increases the cost of its exploitation. These significant drawbacks lead to 
the study of the control of solar photovoltaic water pumping system (SPVWPS), which is the most important 
and appropriate use of photovoltaic in this paper. PV generator, boost converter, DC-link voltage, three phase 
voltage source inverter (VSI), squirrel cage induction motor, centrifugal pump, controller including MPPT, 
rotor flux orientation, rotor speed control unit, and storage tank are the main components of SPVWPS. 

Energy from photovoltaic sources is intermittent and dependent on climatic conditions. Since these 
conditions are not favorable to provide enough energy continuously all day and during the night, 
the photovoltaic generators are generally coupled to the motor-pump unit through a conversion system which 
makes it possible to control the system by optimizing the energy produced by the generator. Therefore this 
energy is exploited to store water in a tank to ensure its availability all times. In addition, the generator has 
non-linear characteristics with maximum power points. These characteristics depend on the irradiance and 
the temperature level. Furthermore, depending on the characteristics of the load supplied by this generator, it 
is possible to find a very large difference between the potential power of the generator and that actually 
transferred to the load in direct connection mode. The first objective of this work is to study maximum power 
point tracking and thus, to transfer of the maximum power of the generator towards the motor-pump drive 
system by means of a boost converter. The second objective is to regulate the rotor flux level at its rated 
value and the mechanical speed of the induction motor at a reference deduced from the maximum power of 
the generator, to extract maximum water quantity, by controlling the VSI. 

In recent decades, many methods for extracting MPP have been developed such as Perturb and 
Observe (P&O) [3-5], Incremental Inductance [3],[6], Fuzzy Logic [7-8], Artificial Neural Network [9], [10] 
Backstepping [11], Backstepping with Integral Action [12], First Order Sliding Mode [13]. 
These techniques differ in many aspects such as sensors requirements, complexity, cost, efficiency, speed 
of convergence, robustness of tracking during irradiation and temperature changes, and necessary 
equipments for digital implementation. P&O has a simple control structure and reduced number of 
measured parameters, however the major drawbacks of P&O algorithm are presence of oscillations around 
the MPP in steady state operation and occasional deviation from the maximum operating point in case of 
rapidly changing atmospheric conditions, such as broken clouds [14]. The control of the induction motor 
remain complex due to the nonlinearity of its mathematical model, high order degree, parameters variations 
and detuning between rotor flux and electromagnetic torque. 

For these reasons, using IRFOC based on classical PI speed controller to implement a high 
performance induction motor drive is often characterized by an overshoot during startup, undershoot while 
applying loads, steady-state error on speed, oscillation in the torque response and overshoot again at load 
removal [15]. Induction motor vector control based on nonlinear approaches, such as Backstepping, 
first order sliding mode, guarantee the global asymptotic stability of the closed-loop system. In addition, 
these nonlinear controllers can solve problems related to classical MPPT. However, these techniques suffer 
from sensitivity to parameters variation, need for measuring all state variables and require observers if 
some state variables are not measurable. To achieve good performance, expensive mechanical sensors 
required in nonlinear control laws and used to measure load torque and rotor speed, are replaced using 
estimators. This paper is organized into five sections: mathematical model of the solar photovoltaic water 
pumping system are described in section II. Section III presents the proposed linear and nonlinear control 
schemes of the system and the developments of estimators for rotor speed, rotor flux and load 
torque disturbances. Section V analyses and discusses simulation results of this comparative study. 
Conclusions and recommendations are given in section VI. 


2. SPVWPS MATHEMATICAL MODEL 

Most of SPVWPS have the added advantage of storing water for use when the sun is not shining, 
eliminating the need for batteries, enhancing simplicity and reducing overall system costs [16]. Figure 1 
illustrates solar-powered water pumping system configuration. Photovoltaic pumping system is composed of: 
— A generator consisting of photovoltaic modules interconnected to form a unit of direct current production. 


— A power conditioning unit, consisting of static converters (boost and inverter), capable of varying 
the frequency and the output voltage according to the available power of the solar generator. 

— Animmersed electric pump unit consisting of an induction electric motor and a centrifugal pump. 

— A hydraulic infrastructure that drives water from its source to a storage tank. 
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Figure 1. Solar-powered water pumping system control scheme 


2.1. PV cell and PV generator mathematical model 

Photovoltaic cell model necessarily involves a judicious choice of equivalent electrical circuits. 
Several mathematical models were developed to represent the nonlinear behaviour of the cell [16]. 
These models differ from each other by the mathematical procedures and the number of parameters. In this 
work, we have chosen PV solar-cell single diode model which is a current source in parallel with one diode 
and both series and parallel resistors. The single-diode model has less computation time and less number of 
unknown parameters [17]. Equivalent electrical circuit is reported in Figure 2. 


Veet 





Figure 2. Equivalent electrical circuit 


The following equations describes the mathematical model of the PV solar-cell 
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Where Zeu» Veen , I 


the solar cell output voltage, photo-current, reverses saturation current, Cell short-circuits current at 298°K 
and 1kW/m2, reference solar irradiation in W/m’, solar irradiation in W/m’, cell’s reference temperature, 
temperature on absolute scale, Cell’s short-circuits current temperature coefficient, electron charge, 
Boltzmann’s constant, and intrinsic series and parallel resistors of the cell respectively. PV generator consists 
of modules connected in series and in parallel to obtain the desired power. Each module is composed of PV 
cells. For the modeling of the PV field, we will therefore start from the basic element which is the PV cell. 
The equations, described PV generator model, are given as: 


one Is, Is , Eref, E, Tn, T, Ki, q, K, Rss and Rsn are the current across the cell, 


Loy = N,I Vpv = N Veen 


cell > 


(4) 
Where and V, are the PV array power and the PV array voltage 


With N, being parallel combinations sets of PV cell and each set consists of N, series-connected 
panels. The SM55 PV panel will be used later in the validation of simulation results. The electrical 
specifications of SM55 PV panel are recapitulated in Table 1. In order to analyze the behavior of the chosen 
photovoltaic panel and PV generator, based on the three-parameter model, we carried out a simulation in 
Matlab/Simulink environment. The results obtained are presented in Figures 3 and 4. Simulation is provided 
to obtain the I-V and P-V time histories of the solar module and PV generator for different solar irradiations 
and temperatures. 


Table 1. Electrical specifications of SM55 PV panel 








Electrical specifications Rating Values 
Maximum power rating Pemax (W) 55 
Rated current at the MPP (A) 3.15 
Rated voltage at the MPP (V) 17.4 
Short circuit current Is (A) 3.45 
Open circuit voltage Voc (V) 21.7 
Current temperature coefficient (mA/°C) 12 

Ns 36 

Np 1 
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Figure 3.Ppv-Vpv Characteristics of PV generator Figure 4. Ipv-Vpv Characteristics of PV generator 


According to simulation results, it is possible to raise the following remarks: When the irriadiation 
increases the voltage and the current increases thus the power increases. In addition, when the temperature 
increases the current is relatively constant, however, the voltage decreases thus power decreases. 
Maximum power of solar module and generator depends on the value of solar irradiation and temperature. 
Components of MPP are summarized in the Table 2 for different values of temperature and solar irradiation. 
These coordinates of MPP will be used to validate the maximum power point tracker. 
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Table 2. Coordinates of maximum power point (MPP) 








Coordinates of MPP TCK) E(W/m) Vinpp (V) Prpp CW) Tryp (A) 
MPP1 298.15 1000 395.80 5007.00 12.65 
MPP2 298.15 500 372.80 2346.00 6.29 
MPP3 323.15 500 326.90 2032.00 6.21 





2.2. Model of induction motor 

Induction motor study consists of looking for all equations connecting input variables to output 
quantities: voltages, stator currents and torque. The different approaches for modelling IM are based on 
electromagnetism and mechanics equations resolution. The model of induction motor remains complex due 
to its high order degree, its strongest nonlinearity and the coupling between flux and torque. To simplify their 
electrical equations, from a six order system to four order system, a dq transformation is used. 
This transformation ensures instantaneous power conservation. Induction motor state space representation, 
expressed in synchrounsly refrence frame, is given by the following [12],[18, 19]. 
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„Tts and R, are dq axis stator current components, dq axis stator voltage components, input control 
components, dq axis rotor flux components, rotor flux, stator electric angular pulsation, rotor electric angular 
pulsation, mechanical speed, stator resistance, stator inductance, electromagnetic torque, load torque, mutual 
inductance, viscous friction coefficient, DC link voltage, leakage coefficient, number of pole pairs, rotor 
time-constant, stator time-constant and rotor resistance. 

This complete model can be considerably simplified by considering an IM supplied by a source 
voltage inverter which is controlled by current supply and by using the indirect rotor field oriented control. 
Indeed, when the stator currents are imposed, the statorique phases dynamics can be neglected. The dynamic 
behavior of the machine is exclusively described by that of the rotor phases. Therefore, in this case, 
simplified second- order model of a current-fed induction motor used to desing simple control laws is 
formulated as [18]: 


d | © Ta, = h0, y Oa, | oa 


dt| ¥ ay, a 0 |Z 


r 


©) 


2.3. Boost converter model 

The DC-DC boost converter is an interface that allows adaptation between the PV generator and 
the other components of the system in order to extract maximum power. It aims to step up the input voltage 
and provides at output an ajdjustable DC voltage. This device contains a diode, a switch, an inductance which 
is added as energy storage element and a capacitor added to output to remove voltage ripple. The following 
state equations describe the boost converter model: 
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dv,, 1-1, 

d G; 

dl, _V,,-(A-w)V, 
dt L, 
dV, (l-u), -I, 
dt C, 


(10) 


2.4. Voltage source inverter model 

DC to AC converter is used to control rotor speed and rotor flux of the induction motor, 
mechanically coupled to centrifugal pump, based on the actual situation of the system. Morover, the inverter 
is able to adjust operating point to an optimal rotor speed deduced from the maximum power 
point. Assuming that the electronic devices of the inverter are ideal and the load fed by the inverter is a 
balanced three-phase load, it is possible to formulate the nonlinear model of this converter based only on 
the switching functions and the DC voltage input. 


Vw] [2 -1 -ifK, 
Vẹ |=2|-1 2 -1 K, V, 
Ven i BA, ie 


(1) 


With K,,K,and K; being the switching functions and V, the voltage output of the boost converter. 
Output currents satisfy the following equation: 


Kiia T Kais E Kiis = 0 (12) 


With i „> andi „being the stator currents of induction motor. 


sa? 


2.5. Centrifugal pump 
Centrifugal pump load torque and mechanical power can be modeled by the folloying expressions [20]: 


ee) 

C, = Ko, (13) 
= 3 

fy EO, (14) 


Where C, is the load torque and P, is the mechanical power. 


3. DESIGN LINEAR AND NONLINEAR CONTROLLERS 

In this section, maximum power point tracking and inverter control will be presented using P&O 
algorithm to control the boost converter. In addition, IRFOC of induction motor based on classical PI speed 
regulator is considered, and the overall system control by two nonlinear control techniques will be discussed. 
Finally, load torque and rotor speed estimators will be developed. 


3.1. P&O algorithm based MPPT technique 

MPPT technique tracks MPP from the measured currents, voltages, or powers. It can react to 
unpredictable changes in PV generator operation. For this, the voltage of the operating point is incremented 
in regular intervals. If the output power is larger, then the direction is maintained for the next step. Otherwise 
it will be reversed. The actual operating point then oscillates around the MPP. P&O disturbs voltage of PV 
generator with AV and observes its PV power variation at the output [21-24]. Figures 5 and 6 illustrate 
flowchart of P&O algorithm based MPPT and its block diagram model in Matlab/Simulink. 
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Figure 5. Flowchart of P&O algorithm based MPPT Figure 6. MPPT block diagram in Matlab/Simulink 


— If AP > 0 then the disturbance of the voltage moves the operating point to a point closer to the MPP and 
the technique continues to disturb the voltage in the same direction. It will move the operating point until 
reaching the MPP. 

— If AP <0 the operating point moves away from the MPP then the voltage is disturbed with an algebraic 
sign opposite to the preceding sign to move the operating point until the MPP is reached. 


3.2. IRFOC based on PI regulator 
Starting with rotor speed and rotor flux tracking errors, we compute first the control laws J e and [ E m 


The second step is to transform these signals expressed in synchrounous reference frame to stationary phase 
coordinate system using dq-to-abc transformation [12], [25-28]. The figure 7 illustrates the block diagram of 
IRFOC based on PI speed. Stator currents references are given by: 


La =Y, pM 
sd rep | (15) 


tek oa. )+K,[(@, ~0,,,)dt 
q P of of (16) 





Figure 7. Block diagram of IRFOC based on PI speed regulator 


4. NONLINEAR CONTROLLERS 
Nonlinear controllers aim to ensure asymptotic convergence of the three outputs to be controlled. 
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4.1. First order sliding-mode control (SMC) 
The real control law u (the duty cycle of the boost converter) is formulated as [29]: 


V-V, Ln 
usd- r etie matj +A f+ 


C, 





sb- msign(s,)— k,s 


E C, L, 


4 (18) 


According to the derivative of Lyapunov function V,, the closed loop system is asymptotically stable: 


V, = $15, = =n |sıl -ks 





(19) 
Define the sliding surface 
ə i=l 
8y(€9,63)= Sa (£) (2 +22] E=E 
: (20) 
The sliding mode inverter control is expressed by: 
=i dæ 
* x 0 Y, = A =y = ref 
P pJ = | a | 027030, dt — msign(sy)— ky 
a, 0 
ag. 
(21) 


From time derivative of second Lyapunov function, it can be seen that the closed loop system is 
asymptotically stable: 


: . 5 
Vy = 5959 = =m |s2l —k,85 (22) 


4.2. Backstepping control 
Let: 
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The MPPT control law u is expressed as[12]: 








L, fE sa Vy 

Shi aa +l 

VG, V, g 
The inverter control I;, and Ka can be given as [30]: 
Tyq = (C43 —4g¥,)/ a, (24) 
* dO, of 

Isg = (C382 + A3 + A30, + )/a¥, 

dt (25) 


It is possible to demonstrate that the time derivative Lyapunouv functions are negative. Consequently, 
asymptotic stability of closed loop system is guaranteed. 


4.3. Rotor speed and load torque estimators 

Load torque disturbance and rotor speed are needed to implement the inverter control. Estimators 
are widely used to avoid expensive mechanical sensors and increase reliability of the overall system. 
These estimators are given by: 


A 2 





C,. = ao, (26) 
aah = aP 1.4 — ay — 430, 
dt (27) 


5. SIMULATION RESULTS 

We tested the linear and nonlinear controllers by simulations using Matlab/Simulink environnement. 
The parameters of the suggested controllers are chosen to achieve high performance. Their sizing is based on 
successive test method. The design parameters of the proposed controllers are chosen as: 

— PlIcontroller: (kK; =5, K,= 200) ; 
— Sliding-mode: (77, =k, = 9000, A, =1200), (77, =k» = 3000) ; 
— Backstepping: (c, =1000,c, =100000,c3 =c, =3000) . 

Solar irradiation variations are chosen in two extreme cases: low irradiation of 500W/m? and high 
irradiation of 1000W/m’. Temperature is fixed at T=298.5°K and is increased to 323.15°K at time t=3.5s. 
maximum power and voltage, rotor speed, rotor flux, load torque and switching functions are reported in 
Figures 8 to 13. As the simulation results show, all proposed strategies achieve the desired objectives. Linear 
and nonlinear control based MPPT algorithm reach maximum power and voltage delivered by the PV 
generator recapitulated in Table 2. IRFOC based on conventional PI speed controller, Backstepping and 
sliding mode force the estimated rotor speed to track the reference value calculated from the obtained MPP, 
and the rotor flux is kept to its desired level. Moreover, these controllers achieve decoupling between rotor 
flux and electromagnetic torque. 

Even if, in transient state, all the controllers present a delay, they respond quickly in steady state. 
Sliding mode and Backstepping show the same behavior: they exhibit fast dynamic performance, stable 
response and have fast and better transient responses than P&O based MPPT algorithm combined with 
IRFOC based on PI speed controller. It can be seen also that the linear controller present oscillatory behavior 
around the MPP in steady state operation, around rotor speed reference and occasional deviation from 
the maximum operating point in case of rapidly changing atmospheric conditions especially in the third 
operating point. However, P&O — IRFOC based on PI speed controller have a simple structure, require only 
measurement of PV voltage and PV current to track the MPP and is independent on the PV generator 
characteristics in the MPPT technique. This technique has the advantage of simplicity that makes it the most 
commonly used in commercial PV system. 
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Figure 12. PV current Figure 13. Load torque estimator 


6. CONCLUSION 

This paper introduced three control strategies for solar photovoltaic water pumping system powered 
by an induction motor driving a centrifugal pump. It provided a comparative study between three control 
strategies of a SPVWPS namely: P&O — IRFOC based on PI speed controller, Backstepping and Sliding- 
mode controls combined with IRFOC. Mathematical models of the main elements of the system were 
proposed. All the controllers reached the MPP under abrupt changing in climatic conditions and tracked 
accurately the desired value of rotor speed and rotor flux. Rotor speed was estimated using only the 
measurement of stator current, the estimate of rotor flux, mathematical model of motion equation and simple 
model of the centrifugal pump. An open loop estimator was proposed to estimate rotor flux. Nonlinear 
controller showed clearly good behavior in transient and steady state, presented fast response time and good 
robustness against variations of climatic conditions and less oscillation around the maximum point in 
the three operating points. Stability of nonlinear controller was guaranteed by means of Lyapunov analysis. 
In contrast, it is not possible to demonstrate stability of P&O — IRFOC based on PI speed controller. In the 
next work, it is possible to propose a robust non linear controller of SPVWPS using a new model of 
centrifugal pump replacing the classical parabolic model. 
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